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Abstract 

In the first part of tliis work (jarXiv: 121 1.65301) it was sliown tliat 
the Hamiltonian formalism of bigravity with the ultralocal potential 
of a general form, under some conditions imposed on the potential, 
had 4 first class constraints and their algebra was the algebra of hy- 
persurface deformations. In this article we suppose that the potential 
also satisfies a special degeneracy condition suggested by de Rham, 
Gabadadze, Tolley (dRGT), and derive the constraint algebra of the 
bigravity and the massive gravity. In contrast to most works in this 
field we do not use neither the explicit form of the dRGT potential nor 
the Hassan- Rosen transformation of variables, but follow some kind 
of axiomatic approach. Also we prefer Kuchaf's Hamiltonian formu- 
lation which involve embedding variables to avoid any noncovariance. 
For the bigravity we get 4 second class constraints, as one degree of 
freedom is included into the phase space in addition to 12 standard 
ones. Given the presence of 4 first class constraints, total number of 
the gravitational degrees of freedom in bigravity occurs 7, and in the 
massive gravity it is 5. 
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1 Introduction 



Despite the fact that the General Relativity (GR) for a very long time stands 
as a Standard Model of gravity, a lot of efforts is undertaken to provide its 
generalizations, and the bigravity theory is one of them. N. Rosen [1] was 
probably the first who coined the term "bimetric" and who advocated a 
theory of gravitation with two spacetime metrics, one of them was a non- 
dynamical, i.e. a background one. Later dynamics for the second metric 
was introduced p], at first such a theory was called "f-g gravity" or "strong 
gravity". The terms "bigravity" and "multigravity" appeared in works by 
J. Kogan and his collaborators [3J where this theory was motivated by the 
boom of extra-dimentional models from the theory side and by the newly 
found effect of acceleration in the Universe expansion from the side of exper- 
iment. It was shown that in bigravity there were two gravitational fields, one 
of them was massless and the other was massive. The massive gravity is a 
long standing challenge for theoreticians. Here S. Deser stays an outstanding 
and permanent sceptic [1]. In a recent article Deser and Waldron [S] claim 
that even if a theory is ghostless, then it is acasual. Also we were reminded [5] 
that this field of research had a rather long history, and remarkably, even the 
most fashionable variant of massive gravity [6] was discussed more than 40 
years ago [7]. Nevertheless we hope that there is some space for the bigrav- 
ity. Both the bigravity, and massive gravity seem to be able to explain the 
accelerated expansion of the Universe, and this is discussed in many papers, 
see, for example, Ref. [8]. 

We think that the de Rham-Gabadadze-ToUey (dRGT) model [6] of mas- 
sive gravity and Hassan- Rosen model of bigravity ^ are very interesting and 
deserve a thorough study. In particular, it is necessary to obtain a more trans- 
parent analysis of its canonical structure, including finding all the constraints, 
calculation of their algebra, classification of these constraints and fixing the 
number of degrees of freedom. Obviously a lot of articles with this work have 
already been published, we mention only a few |9l [TOl IHl [121 IISI EH US], but 
the problem is difficult due to the complexity and bulkiness of the proposed 
derivations coming from the nontrivial form of the potential with the matrix 
square root. 

Let us point out three special features of our method. First, we intend 
to use the benefits of covariant Hamiltonian approach by Kuchar [16], which 
takes into account, in contrast to the famous ADM approach [T7j, not only 
metrics, but also variables describing embedding of a state hypersurface into 
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the spacetime. This allows to vary spacetime coordinates and metric compo- 
nents in the 3 + 1-basis independently. Let us mention that in bigravity we 
have two such bases. Second, we expect that a careful study of bigravity will 
also lead to better understanding of the massive gravity, i.e. of a theory with 
only one dynamical metric. Third, the plan of this article is to replace the 
cumbersome calculations with the dGRT potential by working with a gen- 
eral potential with some axiomatic requirements imposed on it. There are 
two kinds of matter present in formalism, each one interacts (the coupling 
is minimal) with its own metric, and our main results do not depend on the 
concrete form of the matter Lagrangians. We divided our work in two parts 
calling them as the general case and the special case of bigravity. 

In the first part of this work [IB] we obtained the conditions which po- 
tential of a general form U{u,u\riij,'yij) should satisfy in order to have 4 
first class constraints. The second class constraints were used for a construc- 
tion of Dirac brackets. In these brackets the algebra of first class constraints 
was the celebrated algebra of hypersurface deformations (see ADM [T7j and 
Dirac [H]). In counting a number for gravitational degrees of freedom we 
started from a number 12 for two induced metric components rjij, and 
subtracted 4 as a number of first class constraints 71, TZi, so getting 8 as a 
number of free canonically conjugate pairs of variables. This corresponded 
to one massless and one massive gravitational fields, massive one there con- 
tained a ghost degree of freedom. 

In this article, which is the second part of our work considered as a 
whole, we put new restrictions on the potential function in order to meet 
properties of the dRGT potential. Namely, in article [18] it was supposed 
that matrix of second derivatives of the potential over variables u, was 
nondegenerate. Here we require that this matrix is degenerate and has rank 
3. At the same time the potential is to fulfill conditions, found in Ref. [TS] . 
which are necessary for the existence of first class constraint algebra, and one 
new condition on the potential is found in this work. 

To compare our results with the results of other groups one should have in 
mind some difference in terminology. First, we call the theory of two dynam- 
ical metrics as bigravity, following earlier work P] , while some authors call it 
as bimetric theory [20j . Second, we start with treating the lapse and shift as 
canonical variables, so our primary constraints are written as conditions that 
the momenta conjugate to lapse and shift are equal to zero. Some authors 
start with admitting that lapse and shift are Lagrangian multipliers and so 
do not have conjugate momenta at all. Therefore they call our secondary 
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constraints as primary ones, and our tertiary constraint as a secondary one. 
At last, as we start with a wider set of gravitational canonical variables, so 
we finish with a wider set of constraints. The additional pair of constraints 
include one primary and one quaternary constraint which are absent in most 
other publications. The reason for this is that our additional pair of canoni- 
cal variables is replaced by one variable treated as nondynamical, and so our 
primary constraint is absent, whereas our quaternary one is understood as 
an equation fixing this nondynamical variable. 

Really, there are at least three different algorithms which allow to come 
to the results. The first one is to construct Dirac brackets on the base of the 
set of second class constraints. It is suitable that this set is not necessary to 
be complete. The second method is to use more Lagrangian multipliers, but 
work with the Poisson brackets. And the third way is to solve the second 
class constraints (may be implicitly) and express some variables as functions 
of the reduced set of canonical variables. We pay most attention to the first 
method, i.e. in the main text we work with the Dirac brackets, but one can 
see in short how other methods work in Appendices. 

In contrast to the general case studied in article [IB] , as here the potential 
has the declared above new properties, one of the secondary constraints S, 
as well as one of the primary constraints 7r„, change their role, and drop out 
of the primary set of second class constraints (the secondary set arises at the 
end). Therefore, only 6 constraints vr„i,iSi are here involved in the construc- 
tion of the nontrivial Dirac brackets. Then we stand with 6 other constraints 
TTu, S, 71, TZi, and are to study their algebra in the Dirac brackets. The 
calculations show that two brackets {S,7l}i), {S,S}i) look as nonzero on 
the surface of the given constraints. For the first bracket this nonzero result 
leads also to the nonzero bracket S, H^. In order to preserve constraint S in 
the process of Hamiltonian evolution we should treat this result as appear- 
ance of a new constraint Q. This is the celebrated "secondary constraint" of 
Refs. P, [13]. In its turn this constraint has, at first sight, a nonzero Dirac 
bracket with the primary constraint 7r„ present in the Hamiltonian (and so 
depends on unwanted variable u), whereas this bracket remarkably weakly co- 
incides with {iS, S}i). Let us remind that at first, in article by Kluson [12], it 
had been claimed that the last bracket is nonzero, but later it was proved [13] 
that this bracket was zero for the dRGT potential. Here we prove this zero, 
and this is crucial in order to get 7 gravitational degrees of freedom for the 
bigravity and 5 for the massive gravity. We acknowledge that this result, 
missing in the first version of this work [21], has been inspired by article [22]- 
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The content of this article is as follows. 

In Section 2 the notations and results of our previous article [IB] are 
reminded. 

In Section 3 we outline in brief properties of the dRGT potential, inter- 
pret the dRGT conditions in our variables, and discuss consequences of these 
conditions. The 8x8 matrix of Poisson brackets between 8 constraints, con- 
sidered in Ref. [TH| as nondegenerate, now becomes degenerate, consequences 
of this fact are analyzed, and the Dirac brackets are derived for the canonical 
variables. 

The results of Dirac brackets calculations between 6 constraints 7^, T^j, 
7r„, S are given in Section 4. 

In Section 5 the mathematical results on Monge-Ampere equation are 
applied, the crucial identity 6* = is proved, and the Hassan-Rosen trans- 
formation is commented. 

Section 6 contains discussion of quaternary constraint \l/ and Table [2] for 
Dirac brackets between all constraints. The number of gravitational degrees 
of freedom is proved to be 7. 

In Section 7 we discuss the massive gravity, i.e. a theory with only one 
dynamical metric, whereas the second one is a given solution of GR equa- 
tions. In that case there are no constraints TZ, TZi and the Hamiltonian is 
nonzero. Primary constraint is to be included into the Hamiltonian with 
a corresponding Lagrangian multiplier. Secondary constraint S is not neces- 
sary to be included as it appears automatically to guarantee a preservation of 
the primary one in evolution. In its turn to preserve the secondary constraint 
it is necessary to introduce a tertiary one Q. This is a bit different from Q 
which appears in bigravity because it depends on background fields A^, A^* 
and rjij. As in the bigravity case Cl also have zero Dirac bracket with 7r„. The 
Dirac bracket {17,5} is nonzero, so both constraints are second class. The 
quaternary constraint arises as a condition to preserve the tertiary constraint 
in Hamiltonian evolution. The tertiary constraint does not commute with 
the primary one in Dirac brackets, so the Lagrangian multiplier may be fixed 
by the requirement of preservation of the tertiary constraint. This gives 5 
gravitational degrees of freedom for the massive gravity. 

The Conclusion contains a complete list of conditions imposed on the po- 
tential in this article, and our short comments on the related works appeared 
after publication of the first part of this work. 

We prefer in this article to use the same notations as in work [TS]. In 
particular, for indices running from 1 to 3 we use small Greek letters, for 
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indices running from to 3 - small Latin ones, spacetime metrics /^j^ and 

Qf^i, have signature h ++. When the same letter is used for analogous 

quantities constructed with the first or with the second metric, then upper 
bar refers to the second one, i.e. (/^j,. We significantly exploit here results 
obtained in the first article [18] and refer to its formulas in the following way: 
(1. 11) means Eq. (11) of article [TS]. 



2 Kuchaf 's approach in bigravity 



As it is in more detail presented in article [18], the bigravity action is a 
sum of two GR action copies, each given with its own matter as a source of 
gravity, minus a potential 

C = c^f) + C^3)-^U{f,,,g,,). (1) 

^^^^ = T^^v^/'^^^i^ + ^Mir, (2) 
^^'^ = Y^^a'-Rl?^ + c['h<l>\ 9,u). (3) 

where /, g are determinants of the first and second spacetime metrics, R^J} ^ 

i?!f) are their Ricci tensors, G(^) are the gravitational constants, C^S, 

C^lf are Lagrangians of the first and second matter, A are abstract indices 
for these matter fields. The potential 



mf,.,g,.), (4) 

having dimension m^, is constructed as an ultralocal function of tensor 
= 9^^ fua combined of the two spacetime metrics. Then the Hamilto- 
nian of bigravity is also a sum of two copies of GR Hamiltonian plus the 
potential ultralocally depending on components of 3 + 1-decomposition of 
the two metrics. In the Kuchar approach [161 123], where embedding vari- 
ables = e^(r, a;*) are involved in addition to metric, this is done with 
the help of basis (n'^,ef). Here ef = dX^/dx^ are tangential vectors to the 
spatial hypersurface of state, and vector is a unit normal to the hypersur- 
face (future directed). There are two metrics in bigravity, and so we are able 
to construct two future directed unit normals, and respectively two bases. 
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When decomposing a metric in its own basis we get 6 independent compo- 
nents representing the induced metric on the given hypersurface, but when 
decomposing it in another basis we get 4 additional components. We take as 
a default basis constructed with the help of first metric /^j,: 

f^^'^ = -n'^n" + ri'^efe'<, (5) 
g^^ = ^^^^M^- + ^^^(e^^- + eX)+7^^efe,^, (6) 

in another basis denoted as (n'*, ef ), all formulas are changed 

^M- ^ + yigf ej, (7) 

in the same time both induced metrics rjij — /^i^ef 7^- = g^vG-fe'^ stay 
unchanged. As usual, contravariant tensors are inverse matrices f]ijf]^^ = S^, 
lijl''^ — ^i- The matter fields are to be given by their components of 3 + 1- 
decomposition too. Time is given by parameter r, monotonically numerating 
hypersurfaces and not by coordinate in contrast to the ADM notations. 
The components of this time vector field 

Ni" = —- = Nn" + N'ef = Nn" + N'ef, (9) 
dr 

are called lapse and shift. In ADM approach where the embedding variables 
are not used only one coordinate system is at work, and so these functions 
are given as components of spacetime metric in the coordinate basis 

^-7=^' ^^ = ^0. ^=7^' ^« = /o^- (10) 
In order to transform from Kuchaf 's notations to ADM we are to take 

r = x' = X\ (11) 

then 

m = 5i;, ef = 5f . (12) 

Exploiting (n^, ef ) as a default basis to make calculations easier we re- 
place 4 components of the second metric by functions 

1 g-^' 
u^^==, u'^--—, (13) 
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having the simple geometric meaning: u is an inverse of a norm (calculated 
in the second metric) of vector n°, constructed as a unit normal (in the first 
metric) to the hypersurface, and are three projections (calculated in the 
second metric) of coordinate basis vectors onto this unit normal 

u = = (14) 

Then 

N = uN, N' = N' + u'N, (15) 

and in general case studied in Ref. [18] the bigravity Hamiltonian takes a 
following form 



H = J d^x (nU + N'n^ . (16) 



Here expressions 



7^ = n + uU + u'H^ + U, (17) 
TZi = Hi + Hi, (18) 

are first class constraints satisfying in Dirac brackets the well-known hyper- 
surface deformation algebra 



{n{x),n{y)}D = (r/^^x)7^fc(x) + 7]'\y)nu{y)) 5,{x, y), (19) 

{Ri{x),TZk{y)}D = My)ik{x, y) + 7^fc(x)5^(x, y), (20) 
{ni{x),n{y)}D = n{x)6,i{x, y). (21) 



The Dirac brackets were derived in article |18j by inverting matrix of Poisson 
brackets for 8 second class constraints: 

dU - dU 

7r„ = 0, TT,. =0, H + — = 0, + — = 0. (22) 

The first quartet of them arises as primary constraints, and the second one 
as secondary. The lapse and shift functions, as in GR, are Lagrangian mul- 
tipliers standing before the first class constraints. In this work, as in [T8] . 



"'^To be pedantic we should first introduce momenta conjugate to N and iV', and then 
the primary constraints saying these momenta are really zero. So, constraints TZ, TZi arise 
as the secondary ones. The removal of these lapse and shift variables with their momenta 
from the list of canonical variables looks as a gauge, see [1.27] 
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we admit that the potential U = y/fiU{u,u\T]ij,'yij) fulfills the following 
conditions 

dU dU dil 

Q> - 2.^.,.|- - + (,« - - § ^ 0. (24) 

which are necessary for constraints 7^ = 0, T^j to be first class (see Appendix 
A). 



3 dRGT-like potential in bigravity 

To simplify formulas of this article we introduce some new notations: 
-9277 ^ fffj ff.fj 

s - — -n + — s - — -n +— (26) 

K,,{x,y) = {S,{x),S,{y)}, K = L-'KL-\ (27) 

It is claimed [9J that secondary constraints S, Si cannot be solved for lapse 
and shift in massive gravity or in bigravity with dRGT potential. In our 
approach instead of lapse and shift N, iV* we use variables u, u\ so this 
claim now sounds as 



0. (28) 



It is evident from constraints structure fl2Bl) that this Jacobian is really a 
Hessian 

where we give a notation introduced in [18] (1.94). With the new notations 
(!25|) it is possible to write matrix L as follows 



Ui Lij{x) 



(30) 
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Then the degeneracy of this matrix means that for potential of dRGT-hke 
form, which we call in this article as a potential of special form, the following 
condition should be valid 



U, = tj", (31) 



or 

UiU' = U". (32) 

The straightforward check of degeneracy of L and of nondegeneracy of 1 1 Lij \ \ 
for the dRGT potential is problematic. Our method is to suppose that they 
are valid and to consider their consequences. Let us remind that there are 
two other presumptions adopted in this article: 1) that dRGT potential can 
be written as a function of 3 + 1-decomposition components of two metric 
tensors and 2) that this potential fulfills conditions fl23l) . flM|) . which are 
necessary for existence of the first class constraint algebra of hypersurface 
deformations. 

Introduced in our first article [18] matrix L of Poisson brackets between 
8 constraints of the set xa given (l32l) becomes degenerate and as a result 
cannot be used for a construction of Dirac brackets [19j. Two constraints 
from this set will not be included into the first set of second class constraints 
which will be denoted now as xa, A = 1, ..,6: xa, A = 1, ..,6 

Xa = {T^u^^Si) , (33) 

Then the matrix of Poisson brackets for constraints xa will be 

llfcW.x.fe)}ll^(£,,4^, "S")- 

where 

L(x) = \\U,ix)\\, (35) 

K{x,y) = \\K,,{x,y)\\. (36) 

As we suppose that matrix L is invertible than matrix IKxajXb}!! is 
invertible too, and its inverse has a following form 

C(3;,^) = ( i^'\x)Kix,y)l-\y) l-\x)6ix,y) \ _ ^3^^ 
\ — L {x)6{x,y) / 
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Dirac brackets are now defined as follows 

{F,G}^ = {F,G} - jdxj dy{F,XA{x)}C^''{x,y){xB{y),G}, (38) 

or, in more detail form, 

= {F,G}- j dx j dy{F,^:^.{x)}K'^{x,y){7iAy).G}- 
- J dx{F,7ru^{x)}{L-y^{x){Sj{x),G} + 
+ Jdx{F,S,{x)}{L-J^{x){7r^,{x),G}. (39) 

In case both functionals F, G do not depend on u^, these brackets coincide 
with Poisson ones. As second class constraints 

Si = (40) 

can be explicitly solved for u\ and we can work with these variables as with 
functions of variables u, rjij, jij and of Tii, see Appendix [Bl Functions N, 
are Lagrangian multipliers for constraints TZ = 0, TZi = and have zero Dirac 
brackets with all functionals. Primary constraint 7r„ = is to be included 
into the Hamiltonian with a corresponding Lagrangian multiplier 

H = J d^x (nTZ + N'TZi + A7r„) , (41) 

whereas secondary constraint 5 = reappears from the requirement of the 
primary constraint conservation: 

7r„ = {7r„,H}^ = -5 = 0. (42) 

Here we stand with 13 pairs of canonically conjugate gravitational variables 
(m, 7r„), (ryjj, H*-'), (7ij, vr*-'), and with some pairs of matter variables, {iI)a-, H^), 
{(j)AiT^^)-i of course. There is a set of constraints which are weakly zero (i.e. 
they could not be treated as zero before calculation of Dirac brackets): 

7r„ ^ 0, 5^0, 7^ ^ 0, 7^^ ^ 0. (43) 

Next step is to clear whether this set of constraints is full and to separate them 
into first and second class sets. To solve this problem we should calculate 
Dirac brackets between the constraints. 
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Figure 1: Table of Dirac brackets for main variables 



In more detail some Dirac brackets are given below: 

B-„ = -(L-J'(x){nj(x},j„Mh (46) 

4 Dirac brackets for secondary constraints 

In Appendix A we display Poisson brackets between constraints TZ and TZi. 
Replacing these Poisson brackets with Dirac ones we can treat constraints 

7r„i = 0, Si^ 0, (47) 

as strongly equal to zero, i.e. we can omit them and obtain the following 

{7^(x), n{y)}fj = [{rj^'TZk + uu'S)ix) + {rj^'TZk + uu'S)iy)] 5,,{x, y), (48) 

{ni{x),nk{y)}fy = n,{y)6,k{x,y)+nk{x)S4x,y), (49) 
{ni{x),n{y)}f) = n{x)5,i{x, y) + UiS5{x, y). (50) 
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Dirac bracket of constraint S with itself is, for the general potential, not 
zero: 

{S{x),S{y)}f, = {S{x),S{y)} + U\x){S,{x),Sj{y)}W{y)- 
- {S{x),S,{y)}U\y) + {Siy),S,ix)}U\x) = 
= (6^ - WS){x)6,{x,y) - (9^ - U'S){y)6,{y,x), (51) 

where 

^■^(^/'^(^l-S-,,.^)-.^^)^/, (52) 

and a new notation has been introduced 

d d 

Evidently, if the following condition is valid for some quantity Z, 

bZ = 0, (54) 

then this quantity Z does not depend on u when constraints Si = are 
satisfied. 

Dirac bracket for constraints S and TZ is, in general, not weakly equal to 
zero: 

{7^(x), 5(1/)}^ = {u' - uW) S{x)6,{x, y) - {u{U'S), + 5(x, y), 

and so, provides us with a new (tertiary) constraint i7 = which is called 
secondary in Refs. P DS] 

n = [ub-i)[u,n] + b[u,n] + 

- - dU 
+ ({Uk\i + Uiik) + {uk\i + uiik)b)^ + 

+ «.f/-f -(^«'£',■.-['H+yVO|J- 
here square brackets denote coefficients standing before 5-function in the 
related ultralocal Poisson bracket: 

{U{x),'H{y)} = [U,'H]{x)6{x,y), {U{x),'H{y)} = [U,'H]{x)6ix,y), 
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and indices of quantities m*, W are moved down by means of metric 7jj. First 
class constraints TZi as usual generate spatial diffeomorphisms 

{5(x),7^i(?/)}^ = -S{y)6,^{y,x), 

{Q{x),n,{y)}^ = ~Q{y)6,{y,x) (56) 

It is possible also to calculate the following brackets 

{n{x),nuiy)}i) = S{x)6{x,y), 

{fi(x),7r„(t/)}^ = e\x)5,ix,y)-&{y)5,iy,x). (57) 

In Ref. |13] it was proved that S commutes with itself for the dRGT po- 
tential. So, it is reasonable to hope that our potential should satisfy condition 
= 0, i.e. 

(u^o(si-2,,,^y,''l-^U = (58) 

The two requirements might be fulfilled at once by this: 

1. the tertiary (often called secondary) constraint will appear; 

2. this appearing constraint will not depend on variable u. 

Eq. (158|) will be proved in the next Section. It is evident that if B* were 
nonzero, then = would be an equation fixing variable u and there would 
be 3 second class and 4 first class constraints for 13 pairs of gravitational 
canonical variables in bigravity. This would give 7+1/2 gravitational degrees 
of freedom for the bigravity, and 5 + 1/2 for the massive gravity in line with 
the conclusion first made by Kluson [12]. 



5 Transformation of variables 

Previous sections were written before we discovered Comelli et al articles [22| 
[2l] . Their analysis refers to the case of the massive gravity, the bigravity is 
not considered. Here we try to apply their ideas to the bigravity case. The 
main point is to use the results on Monge-Ampere equation obtained by 
mathematicians [25]. First, we can see that (1, — [/*) is a null- vector of the 
Hessian matrix 
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so, in order to be close to article [2S], we introduce a new notation: 
It is suitable to replace potential U by function ^('C\%)7ij)) such that 



dil dV 



dU 



then it follows that 



Integrability condition for Monge- Ampere equation 



(59) 



(60) 



Det 



du'^du^ 







(61) 



is as follows (see [25] for details) 



(62) 



At last, it is proved in Ref. |25] that any implicit solution for system of 
equations ( 16T1) can be written as follows [25] 



(63) 



1 udW 

where W = Vl^(^*, 7ij) is an arbitrary function, and 



ik 



If we apply differential operator D to Eq. (B^ . and take into account 
and fl60|) we get the following 



dV 

- 

(64) 
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Now 



then after repeating this operation and dividing the result by factor 2 we 
obtain 

2e7.^^l^ + ^'l.^iDf^ - eV - ^''^ = 0. (65) 

d-fke d-fke d-fke ' 
as this expression depends on -u, -u* only through variable ^* and given Eq. ( l62ll . 
it is annihilated by operator D, so Eq. ( 158|) is proved. 

Here it is interesting to compare our results with the argumentation given 
in the pioneering work [10]. Variables here play the role analogous to A^* 
there, u is similar to A^, and variables ^* or f7* are analogous to n*. Two of 
the three criteria for the existence of the Hamiltonian constraint formulated 
in Ref. [10] are fulfilled here as one can see from Eqs. f l59|) and fl63|) . It is 
proved in Ref. [10] that their third criterion is fulfilled automatically when 
the other two are satisfied. And let us consider a function 

du d^^ 
on the surface of constraints Si, then 



du 



dV 

v-cn„ — ^0, 



Si=0 

therefore on this constraint surface we have 

~ 0, 

as it has been stated in Ref. 1101. 



6 Dirac brackets for tertiary and quaternary 
constraints 

For the potential satisfying 0* = we are to continue the process of cal- 
culating Dirac brackets between constraints in looking for new constraints. 
Unfortunately, expressions for constraints evidently become more compli- 
cated with each step. We derived tertiary constraint Q above, see Eq. (^^. 
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Figure 2: Table of Dirac brackets for constraints 



In order to preserve this constraint in the process of evolution we are to 
ensure the following condition 

fi = {i],H}^ = I N{n,n}f, + N'{n,n^}f, + x{n,7Tu}n^o. (67) 

As the second and the third Dirac brackets here are zero on the surface of 
already known constraints, due to Eqs. (156|) . (1571) it is necessary to calculate 
only the first bracket: 

{^(a:),7^(|/)}^ = v]>(x)(5(x,l/), (68) 
that will give us new (quaternary) constraint 

^ ^ 0. (69) 

The Jacobi identity 

{n, {7^„,7^}^}^ + {7^, {i],7r„}^}^ + {ttu, {7^,^]}^}^ = o, 

if we take into account Eqs. flFTl) . gives us the following result 

{7r„(y),vl/(x)}^5(x,z) = -{n{x),Siz)}j,5iz,y), (70) 

and so it is evident that if {Q, S} fj is nonzero, then {vr^, ^ is nonzero also. 
But it is easy to check that there are nonzero terms in {Q, 5}^, for example, 

4/^(/)/£(9) / n\ d'^V f Tx\ 

[[V, H], = Ylij - Tjij— Tlmn " Inm- , (71) 



Vl \ 2 / dr]ijd'Jmn 
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that cannot be cancelled by other terms and do not appear in other con- 
straints. Therefore both Dirac brackets in Eq. fl7Up are nonzero. Of course, in 
looking for the so-called partially massless case, see Refs. [261 ISZl |28l |29| 130]. 
one needs a more detailed study. 

There is a potential problem if \E'(x) occurs not a function, but a differ- 
ential operator, even after taking into account all constraints. Then Eq. (1571) 
can be solved for u, though the solution will depend on A^. 

We combine schematically all the results of our calculations of the Dirac 
brackets in Table [21 these results are displayed above in this text in more 
detail. One can see that there are 3 blocks standing on the diagonal: two 
nondegenerate 2x2 matices for two pairs of the second class constraints 
(vr^, and {^,S) and one degenerate 4x4 matrix for the first class con- 
straints {TZ, TZi). Therefore, for 13 pairs of the gravitational variables we have 
4 first class constraints and 4 second class constraints, then this corresponds 
to 7 degrees of freedom for the bigravity. 



7 Hamiltonian approach to massive gravity 

Let us consider now the case when the first metric is a fixed GR solution 
(maybe, with a cosmological constant) and only the second one is dynami- 
cal. We call such theories as the massive gravity or the bimetrical gravity 
following N. Rosen [I]. Recently they have got a revival [211 E21 ESI El]- 
Given metric /^^ and embedding variables e°(t, x^) we are to treat N{t,x^), 
N^{t,x''), riij{t,x^) as the known functions. For the Hamiltonian we can use 
the following expression 

H = J £x [n (uH + u'Hi + U)+ N'Hi + Att^ + AV„, ) . (72) 

It is similar to bigravity Hamiltonian given that "H, "Hj are now zero, as metric 
f^j_i, is now a fixed solution of GR equations. Of course, the Poisson brackets 
now do not involve variational derivatives over rjij. There are 4 primary 
constraints 

-Ku = 0, 7r„« = 0. (73) 
and, in their turn, they generate 4 secondary constraints 

- dU - dU 

n + — ^S = 0, 'H. + — ^S, = 0, (74) 
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where U = U{u, u\ 'jij, rjij). The Poisson brackets matrix of these constraints 
should be degenerate for dRGT potential as in the bigravity case. Then we 
again construct the Dirac brackets (of course, they are a bit more compact 
than in the bigravity, as the Poisson brackets now do not involve variational 
derivatives over by inverting 6x6 matrix constructed of the second class 
constraints = 0. As the Hamiltonian is a first class quantity we can drop 
out the second class constraints to simplify its form 

Hreduced = / d'x (^N (uH + u) - N'^^ + Att^^ . (75) 

Primary constraint 7r„ is still present in the reduced Hamiltonian accompa- 
nied by its Lagrangian multiplier A, and secondary constraint S reappears 
as a condition of preservation of 7r„ in evolution 

T^u = {tT-u, Hrcduccdj^ = -NS = 0. (76) 

Constraint S can be solved for the ghost variable. For consistency of the 
evolution we should calculate the Dirac bracket of the secondary constraint 
with the Hamiltonian. Applying Eq. ( l38l) we can write this Dirac bracket 
through the Poisson ones 

{S{x), Hreducedjo = {<S{x), Hreduced} " U^{x){Si{x), Hrcduccd} — 

' d'y {u\x){S.{x),Sj{y)} + {S{x),S,{y)}) N'{y). 
The result is a new expression 

{S{x), Hreducedlc = ^, (77) 

which is necessary to be zero 

n = N (uD -i)[u,n] + 

_ _ _ _ afr 



kl 



- f)TI r)TJ 
(^■M + 7V.y'^^ + 2W^7,.g =0. (78) 
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Therefore, we have got a new constraint. Given zero Dirac bracket {Cl, vTu}^, 
this constraint does not depend on variable u when Si constraints are taken 
into account. Eq. flTHj) is hnear in momenta tt*-' and so can be solved for the 
momentum variable which is conjugate to the ghost variable. 

Next step is to derive the Hamiltonian equation for constraint fl 

n = {n, Hreduced}/) = ^, 

and to require that Q should be zero. So we get the last constraint \E', which 
serves for fixing variable u. At the end, we have 4 second class constraints 
S, Q, TTu, \E' for 7 pairs of gravitational canonical variables. Two second 
class constraints ^ and ir^ serve to kill a pair of canonical variables (u,7ru). 
Two other second class constraints serve to kill the ghost degree of freedom. 
Therefore this gives 5 degrees of freedom for massive gravity. 

8 Conclusion 

We have studied here the canonical structure of a model which, as we hope, 
should be isomorphic to widely discussed dRGT theory [6l HI [TOl [TTl [121 [131 
[13] of massive gravity and bigravity. In this work the potential should satisfy 
4 conditions: 

1. the potential can be expressed as a function of 3 + 1-decomposed metric 
variables U{u,u\ rjij, ■jij); 

2. the potential allows 4 first class constraints TZ, TZi in bigravity case; 

3. the Hessian of the potential in the lapse-like and shift-like variables 
M, is degenerat^; 

4. the Hessian in the shift-like variables is non-degenerate. 

We conclude that under these suppositions we get a theory with 5 gravita- 
tional degrees of freedom in the massive gravity case and with 7 degrees of 
freedom in the bigravity case. It would be interesting to prove whether the 
dRGT-potential is a unique realization of the above axioms. 

^Our variables can be expressed through lapses and shifts of the two metrics as follows 
u = N/N, = {N' - N')/N, see Eq.^. 
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Since the first part of this work [IB] has appeared, some interesting and 
related articles have been published. Kluson's preprint [20], in fact, was 
simultaneous to [18] , and the same algebra of first class constraints in bi- 
gravity was proved under some conditions imposed on the potential, the 
dRGT case was not included. Distinctions from our work are in more narrow 
class of potentials, in the choice of variables, and in using ADM approach, 
not Kuchar's one. Alexandrov et al [13] developed Ashtekar-like approach 
and variables. These authors were careful to mention some difficulties with 
the reality conditions. Their main conclusion is a confirmation of statements 
made in articles [9], at least for one variant of the dRGT potential. Many 
interesting ideas are suggested, including a choice of an arithmetic average 
of triads related to different metrics. In article by Deser and Waldron [S] a 
simple derivation of Hamiltonian constraint is demonstrated, the interesting 
reference to considerably earlier work is given, and a statement on superlu- 
minal effects in massive gravity is made. Certainly, to make final conclusions 
on the consistency of any theory with two or more spacetime metrics it is 
necessary to provide a deep analysis of causality problem, but we here do not 
plan to discuss it. Two other articles by Kluson [351 EH] have appeared later 
and contain some pessimistic conclusions related to the bigravity. In contrast 
to our results it is claimed that constraint playing the role TZ becomes second 
class and so the theory loses the diagonal diffeomorphism invariance. As we 
have shown here, TZ stays the first class constraint, as it is in GR. Let us 
mention that calculations are rather long in this theory, and so a choice of 
variables matters. The choice made in Ref. [35] is in line with democracy 
stereotype whereas we give a preference to one metric over the other (the 
revolution is allowed, see (1.59)), this maybe the reason for the discrepancy. 
The appearance of e-print by Comelli at all [22] was important for us. The 
similar plan to study the massive gravity with the potential of a general form 
was proposed there (the bigravity case was not discussed). Reading Ref. [22] 
moved us to revise the conclusions made in the first version of this work [21] . 
Unfortunately in preparing [21] we were not informed on the preceding pub- 
lication by that group [21]. 
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A Calculating Poisson brackets of IZ, IZi 

In calculating Poisson brackets between constraints TZ, TZi we can treat m, 
as functions, and not as canonical variables, because their conjugate momenta 
do not appear in these relations. Then potential U has nonzero Poisson 
brackets only with gravitational momenta H*-' , vr*-' , this results in appearance 
of derivatives dU/drjij and dU/d'jij. In derivation of the first bracket 

{ni{x),n,{y)} = ni{y)6,,{x,y)+n,{x)64x,y), (79) 

we take into account that Hi and "Hj commutes, and each of them satisfies 
Eq. (149|) . then (1791) is valid. By straightforward calculations we get 

{n{x),n{y)} = {p'^Uj + uu'S - {ri'J - - u'u^)Sj + Q*) {x)b^^x, y) - 
- (rf^n^ + uu'S - (rf^ - u^-i'^ - u'v?^ Sj + Q') {y)6^i{y, x) ^ 
^ {v''n, + Q') ix)S,{x, y) - {r^'^n, + Q") {y)5,{y, x), (80) 

and also, 

{7^^(x),7^(?/)} = TZ{x)5^i{x,y) + u^iS5{x,y) + 
d 

+ —(Q{{x)6{x,y)) + 
d 

+ — (u^Si{x)6{x, y)) + ul,Sj6{x.y) ^ 

^ n{x)5,{x,y) + ^(Ql{x)5{x,y)), (81) 

where simpler expressions arise after taking into account constraint equations 
5 ~ and 5i ~ 0. Quantities 

dU dU dU 

Ql^2,,,— ^2,,.— -u'-^-SlV, (82) 

0..,„.,^,||_„-„|,(,,«.„V'-.V)§. (83) 

should be equal to zero in order TZ, TZi to be first class constraints. 

As has been observed in articles [201 [T5] . Eq. ( ISTll becomes simpler if one 
add to TZi a combination of primary constraints: 

1Zi=TZi + -KuU^i + TT^fcM^.i + (7r„«M^),fc, (84) 
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then we obtain 

{%{x),n{y)} = 7^(x)5 y) + ^ {QKx)5ix, y)) , (85) 

but with the replacement of TZi for TZi the primary constraints appear in 
Eq. ([HQ]). 



B Solving Si constraints 

Second class constraints may be used not only for deriving of Dirac brack- 
ets, but also for straightforward exclusion of conjugate variables {u, tTu) and 
(u'^jTiui) from the Hamiltonian. When the chosen potential is simple enough 
(see, for example, [37]), it is possible to express 3 functions (or g^^) through 
the rest variables after solving constraints 

Hi + dU/du' = 0. (86) 

Given these solutions 

u'' = u\u,T]ij,-fij,ni), (87) 

and then substituting them into constraint + dU/du = 0, we can express 
also u (or g^'^) as a function of the rest variables 

u = u{r]ij, 'jij, ni'jij, n'^cpA, TT^))- (88) 

It is stated that for potential of a special form proposed and studied in 
articles [6l |9] it is impossible to solve the secondary constraints and 
express m as a function ( !88l) . As above, we suppose that potential U has such 
a property that 3 constraints Si can be solved for variables u'^, matrix \ \Lij\\ 
is invertible in this case. By substituting these solutions back into equations 
Si = we evidently get identities, and then 



du 

or equivalently. 



^54.^,. ^0, (89) 



dudu^ du^du'' du 
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From this we find 



du" 'L-^fX" . ^^-(L-^y\f" ■ (92) 



dr]mn ^ ' du^dr]mn d-fmn ^ ' Ou^O^., 

Next, after substitution of the derived expressions u'' = u^{u, rjij, 7,^, T-Li) into 
the forth constraint, we get an expression, which according to supposition 
made in articles [6l |9] does not depend on u: 



'-'reduced TL ~r „ 
OU 

and then 



(93) 



a>Sreduced _ d^U ^ de d'U _ ^dU _ ^ 

du du^ du dudu^ du 
In order to commute the constraints we can use equivalently the canonical 
Poisson brackets instead of Dirac brackets, if we take into account dependence 
given by Eq. ([8] 



C Analysis of bigravity with Lagrangian mul- 
tipliers method 

Apart from the two methods used above (Dirac brackets construction in Sec- 
tion 3 and solving of constraints in Appendix [B]) there is an approach based 
on calculation of Lagrangian multipliers. Here we start with the Hamilto- 
nian containing all primary constraints with their Lagrangian multipliers, i.e. 
functions which should to be varied to get constraints as equations of motion, 
but which have zero Poisson brackets with any variable, 

H = y" [nTZ + N'TZi + X-Ku + AV„,) . (95) 

Then we check for consistency of 4 primary constraints 

r dTZ 

TTu = {tTu, H} = y £xN{7ru, 7^} = -— = ~NS = 0, (96) 

r dTZ 

= {n^^,B.}= d'xN{7r^^,n} = -— = -NSi = 0, (97) 
J ow 
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and get 4 secondary constraints 



Si ^ 0. 



(98) 



We should not add secondary constraints to the Hamiltonian, but we are to 
check their evolution equations for consistency 

S = y'iV{5,7^} + iV'^{5,7^4 + A{5,7^4 + A^{5,7^„.} = 0, (99) 
S^ = y'A^{5„7^} + A^'={5„7^fe} + A{5„7^4 + A'={5„7r„.} = 0,(100) 

where omitting 5-function we can write 

{Si, vr^fc} = Lik, {Si, 7r„} = Ui = {S, 7r„»}, {S, 7r„} = JJ". (101) 

As Lik is invertible matrix we can solve Eq. f llOOp for A^ 

A*= = -XU'' - (L-^)''^ (^N{Sj,n} + N'{Sj, TZi}) , (102) 

and then substitute this into Eq. (1^ . Evidently, multiplier A is canceled 
due to Eq. (IHTi) . and therefore Eq. (1^ is a new constraint (compare with 
Eq. dSSD): 

S ^ Nn^O, or n^O. (103) 
Next step is to check a consistency of this last constraint 



N{n, 7^} + N''{n, Uk} + \{n, n^} + \''{n, tt^^} 



N 



{n,n}-{n,n^.}{L-y^{s„n} 



+ N' 



{n,n,} - {Q,n^.}{L-')''^{s„n,} 



+ 



+ \[{n,7r^}-U'{Q,7r^.} = 

= I N{n,n}f, + N'{n,ni}f, + x{n,TTu}D^o. 



(104) 



We can see from Eq. (IS7|) that the last Dirac bracket is proportional to 9*, 
so it is zero. Then Eq. fll04p gives a new constraint. As f2 is a spatial scalar 
we have 

{n,nk}D^o, (105) 
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and therefore the new condition is as follows 

y" A^{^], 7^}^ = ^ ^ 0. (106) 

Next we should check 

^ = {^,H} = y'iV{^,7^} + iV^{$,7^J + A{^,7^„} + A'={$,7r„4 = 
= y'iV{^,7^}^ + iVH*,^fc}£) + A{^,7^J^ = 0. (107) 

As, according to Eqs. ([70]), fITB . 

{vl/,7rj^^0, (108) 
then Eq. fll07p can be solved for A and the analysis of constraints is finished. 
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